EfFective action for superfluid Fermi systems in the strong-coupling limit 
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We derive the low-energy effective action for three-dimensional superfluid Fermi systems in the 
strong-coupling limit, where superfluidity originates from Bose-Einstein condensation of composite 
bosons. Taking into account density and pairing fluctuations on the same footing, we show that 
the effective action involves only the fermion density pr and its conjugate variable, the phase Or 
of the pairing order parameter Ar. We recover the standard action of a Bose superfluid of density 
Pr/2, where the bosons have a mass = 2m and interact via a repulsive contact potential with 
amplitude qb = 4TTaB /rriB , ub = 2a {a the s-wave scattering length associated to the fermion- 
fermion interaction in vacuum). For lattice models, the derivation of the effective action is based on 
the mapping of the attractive Hubbard model onto the Heisenberg model in a uniform magnetic fleld, 
and a coherent state path integral representation of the partition function. The effective description 
of the Fermi superfluid in the strong-coupling limit is a Bose-Hubbard model with an intersite 
hopping amplitude tB = J/2 and an on-site repulsive interaction Ub = 2,Jz, where J — At^ /U {t 
and —U are the intersite hopping amplitude and the on-site attraction in the (fermionic) Hubbard 
model, z the number of nearest-neighbor sites). 

PACS numbers: 74.20.Fg, 71.10.Fd, 05.30.Jp 
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INTRODUCTION 



Recent progress in the experimental control of ultra- 
cold atomic Fermiii2i2iii& gases has revived the interest in 
the crossover from the weak-coupling BCS limit of super- 
fluid fermions to the strong-coupling limit of condensed 
composite bosonsiSii In this paper, we derive the low- 
energy effective action for a superfluid Fermi system in 
the strong-coupling limit, both in continuum and lattice 
models. The latter may be relevant for high-Tc super- 
conductors or ultracold Fermi gases in an optical lattice. 

A Bose superfluid is described by a complex fleld 
'fpBr — y/PBrG^^^' whcrc pBr is the bosou density at po- 
sition r in space. The equation of motion derived from 
the standard action of a Bose system leads to the Gross- 
Pitaevskii equation, i.e. a non-linear Schrodinger 
equation for the ipB fleld. The Gross-Pitaevskii equa- 
tion yields a simple description of quantum macro- 
scopic phenomena Hke the Josephson effect or the flux 



quantization 
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and has proven to be a tool of choice 



for the understanding of many phenomena in ultracold 
atomic Bose gasesjiS In Fermi systems, there is in gen- 
eral no simple relation between the amplitude of the su- 
perfluid (pairing) order parameter Ar and the fermion 
density pr- This suggests that a minimal description, 
aiming at making contact with the standard description 
of a Bose superfluid, should at least include the super- 
fluid order parameter Ar and the density pr from the 
outset. In the strong-coupling limit, where superfluid- 
ity originates from Bose-Einstein condensation (BEG) of 
composite bosons, we expect the description in terms of 
Pr and Ar = |Ar|e*^'' to be redundant and the superfluid 



to be described by a single complex fleld ^Ar = ypTj^c}^" 
{pr/'^ being the density of composite bosons). 

Previous studies of the BCS-BEG crossover in super- 
fluid Fermi systems can be divided into two categories. In 
the flrst type of approach )<^i^i?i^"^i^i^i^^i^T|ifi the density Pr is 
not considered explicitely and a pairing fleld A^^ is intro- 
duced by means of a Hubbard-Stratonovich transforma- 
tion of the fermion-fermion interaction. In the BEG limit, 
the standard action S[ij)* ^ij)] of a Bose superfluid is re- 
covered if one identifles ij^r to A^^ (after a proper rescal- 
ing). For a continuum model, the bosons have a mass 
ijiB — 2m and interact via a repulsive contact potential 
with amplitude qb = ^iraB/mB, aB = 2a (a is the s- 
wave scattering length associated to the fermion-fermion 
interaction in vacuum) . The main (conceptual) difficulty 
of this approach is that the Hubbard-Stratonovich field 
A^^ is not the physical pairing field Ar = |Ar|e*^'' but 
rather its conjugate fieldfi^ Although both fields coin- 
cide at the mean-field level, they differ when fiuctuations 
are taken into account. As a result, V'r oc A^^ does not 
correspond to pr/2e'^"' as expected. 

In the second type of approach to the BGS-BEC 
crossover (iSiSSi21 the physical density and pairing fields, 
Pr and Ar, are introduced from the outset. For con- 
tinuum models, only the weak-coupling limit has been 
considered piSi21 For lattice (Hubbard) models in the 
strong-coupling low-density limit, one finds that the or- 
der parameter am plitud e and the density are tied by the 
relation |Ar| = pr/2, so that the low-energy effective 
action c an be written in terms of a single complex field, 
Ar = •\/pr/2e'^'' = -0ri^^^ In the continuum limit, one 
finds that the (composite) bosons have a mass ras = 1/J 
and interact via a repulsive contact potential with ampli- 
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tude gB = 8J (in two dimensions), where J = At^ /U (t 
being the intersite hopping amplitude and {U > 0) 
the on-site attractive interaction) 31' 

Most of the theoretical works on the BCS-BEC 
crossover in ultracold atomic Fermi gases have been for- 
mulated within a fermion-boson model^ aiming at in- 
corporating the molecular states involved in the Fesh- 
bach resonance which drives the crossover. While the 
equivalence of the fermion-boson model to an effec- 
tive single-channel model in the crossover region may 
be questionabler^iSi both models are equivalent in the 
strong-coupling limit. 

The outline of the paper is as follows. In Sec. lUl we 
extend the approach of Ref. to the strong-coupling 
limit of a continuum model. The particle-particle and 
particle-hole channels are considered on the same foot- 
ing, and the (physical) density {pr) and pairing (Ar) 
fields are introduced from the outset. The low-energy 
effective action is derived by assuming small fluctuations 
of the collective fields about their mean-field values. We 
find that fiuctuations of pr and Ar are not independent, 
so that the low-energy action can b e written in term of 
a single complex field tpr = \J Pr/Se*^"'. We recover the 
standard action of a Bose superfiuid with ms = 2m and 
qb — A-KaB/ruB^ ttB — 2a. For a lattice model fSec. lIIIII . 
we follow the approach introduced in Ref. We map 
the attractive Hubbard model onto the half-filled repul- 
sive Hubbard model in a uniform magnetic field coupled 
to the fermion spins. In the strong-coupling limit, the lat- 
ter reduces to the Heisenberg model in a uniform field. 
The low-energy effective action of the attractive model is 
finally deduced from the coherent state path integral rep- 
resentation of the Heisenberg model. In the low-density 
limit, where the Pauli principle (which prevents two com- 
posite bosons to occupy the same site) should not matter, 
|Ar| ~ \/pr/2 and the superfiuid Fermi system can be 
described by the complex field V'r = \/ Pr/2e*^"' . We find 
that the effective description of the Fermi superfiuid is 
a Bose-Hubbard model with intersite hopping amplitude 
tB = >//2 and an on-site repulsive interaction Ub — ^Jz 
(where z is the number of nearest-neighbor sites). 



II. CONTINUUM MODEL 

We consider a three-dimensional superfiuid fermion 
system with the action = 5*0 4- ^int , 



2m 



•S'int, — 



where ci*J are Grassmann variables, Cr — (crt,Cr|)-^, r 
an imaginary time and /? = 1/T the inverse temperature. 
—g is the attractive interaction between fermions {g > 0). 
The chemical potential p fixes the average fermion den- 
sity po- To suppress ultraviolet divergences appearing 



(1) 



in the perturbation theory, one regularizesi^ the local 
fermion-fermion interaction with a cutoff A acting on the 
fermion dispersion: ek — |kp/2m < A. g and A de- 
termine the s-wave scattering length a defined by the 
low-energy limit of the two-body problem in vacuum. 



m 



ek<A 



d^k 1 

(2^2^- 



(2) 



a is negative for small g and diverges when g — 27r^/mA. 
For g > 27r^/mA, there is a two-body bound-state (com- 
posite boson) with energy Eb = — and the scat- 
tering length a is positive. The latter also determines 
the extension of the bound-state. Low-energy properties 
depend solely on a (and not g or A); we shall therefore 
take the limit g ^ and A ^ oo with a fixed. In the fol- 
lowing, we consider the BEC limit defined by poa^ <C 1 
(a > 0) , where superfiuidity originates from BEC of com- 
posite bosons. 

The (real) density and (complex) pairing fields. 



Pr 
Ar 

A! 



— CrCr, 

= cJcr^Cr, 

= CrXCrT, 

— Cr|Cr|, 



(3) 



can be introduced in the action by means of real 
{p^^,p^^) and complex (A^^) Lagrange multipliers: 

-57|Arp-f (p?-5r^^) 

+ipf^{pr - 4cr) + ipP(S'r - ctfl'^Cr) 

-hz[aP(a;-c;^c:j + c.c.]|. (4) 

{a^,ay,a^) denotes the Pauli matrices. Integrating over 
p^^,p^^, A^^ and pr, S'r , Ar, we recover the original ac- 
tion (QJ if we choose a + j = 1^ The relative weights 
a and 7 of the particle-hole and particle-particle chan- 
nels are arbitrary. All the resulting effective actions are 
equivalent when treated exactly. However, to recover the 
mean-field results from a saddle-point approximation, we 
take a — J — 1. When only low-energy long- wavelength 
fiuctuations about the mean-field state are considered, 
there is no overlapping of the two channels and there- 
fore no overcountingjJ^ Note that by integrating out the 
physical fields S^, pr and Ar, one recovers the action 
S[c, p^^, pj?^, Aj^^] which is generally obtained by means 
of a Hubbard-Stratonovich decoupHng of the interaction 
term. Thus the Lagrange multipliers p^^, p^^ and A^^ 
enforcing the constraints Q can also be seen as Hubbard- 
Stratonovich fieldSflS In the following, we neglect spin 
fiuctuations (p^^ and S^) since they do not play an im- 
portant role when the interaction is attractive. 
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A. Mean-field theory 

The mean-field theory is obtained from a saddle-point 
approximation where the fields Pr, ^r, P^^ and A^^ 
are taken space and time independent. The saddle-point 
equations read 



Ao = (CrXCrt), 
— {C„t C„ 1 ) , 



zAf = .gAo, 



9^0- 



(5) 



With no loss of generality, we can take Aq = Ag real. 
zAg ^ = zAg ^ is then real at the saddle point. It is con- 

A^^* (so 



venient to redefine zAf? ^ 



A^s and iA}}^* 



that Ag ^ = A^^ is real) and absorb ip^^ in the defini- 
tion of the chemical potential. The mean-field action is 
then (up to an additive constant) 



Smf — 



9r - M - ^ I Cr 



-Af'^icl^^i^+cc.) 



(6) 



From we readily obtain the normal and anomalous 
Green functions, 

-iuj - 



G(k, iw) = — (co-(k, zct;)c* (k, zcj)) 



i^(k, iw) = — (co-(k, zci;)cg-(— k, — icj)) 



AgHS 



F*{k,iu;) = -(c;(-k,-zc^)c:(k,ic^)) =i^(k,zc^),(7) 



where = (Ck + ) ' > = Ck - P, and ct = ~a. 
Co- (k, iw) is the Fourier transformed field of Cra and uj a 
fermionic Matsubara frequency. Using ^ and Q, we 
can rewrite the saddle-point equations J^J as 



m 
47ra 



1 



(8) 



1 

'k \2ek 2£^k, 

Av ^k, 

where = J ci^k/(27r)^. Eqs. (jHl determine the chem- 
ical potential p and the order parameter Ag^ = gAg. 
In the strong-coupling limit poa^ <C 1, one obtains (see 
Appendix IB| 

^ r(l-2w'), 



^0 — 



47rpo 



1/2 



1 + -poa' 



(9) 



B. Low-energy effective action 

In this section, we derive the low-energy effective ac- 
tion for the physical fields Pr and Ap. Since our deriva- 
tion partially follows Ref. Ilil we describe only the main 



steps (technical details are given in Appendix EJ • The 
main assumption is that the collective bosonic fields pr, 
p^^ , Ar and A^^ weakly fiuctuate about their mean-field 
values. 

Starting from the action Q (with a = 7 = 1), where 



Ar - |Ar|e'^^ 
we perform the change of variables 

Cr^Cre*^^ A^'^-^A 



(10) 



(11) 



oHS 



Af -Kz 



A?^* (recall that a 



We then consider the shift 
Af^ -f iAf^ and iA^' 

factor i has been included in A[/^ and Ag ^ ) , so that the 
Hubbard-Stratonovich fields and A^^ now describe 
(small) fiuctuations about the mean-field values. This 
leads to the action 



S'mf + 



dT \ d^V 



8m 2 
iAP%zAf -^2AHS)|A 



V6'r • V 

i{Afcl^cl^^c.c) 



Am 

,HS„* 



9 2 



-g\A^\'-■^pi + i^p^ 



„HS 



ipP)pr 



(12) 



where V = V — V. Here we write the chemical potential 
as p = pmf + Mb/2 where pmf is the chemical potential 
in the mean-field approximation. The next step is to shift 
ier/2 + (V6lr)V8m - AiB/2, and to 



pHS, zpHS 



ip, 



,HS 



introduce Nambu spinors (pr = (crtjC*j^) . This gives 



S = 



Smf + S' 

HS 



-Pr(«Po^ 



dT J d^r 

HS „■ A HS 
HS 



-5|Ar 





(2A;^^+zA«^+zA^^'-)|Ar 



(V^r)' MB\ 

8m 2 ) 



(13) 



where 







S' = 












ior ~ 


(j)lT^(j)r = cjc 




0];r+(/)r = c* 


.?0r = 


^Jr^l/ir = Cr 




--^0tV0r 

2m 



, HS* •- 



(14) 



rT'-ri: 



2m 



-cIVCr 



(15) 



{Tx,Ty,Tz) are Pauli matrices acting in Nambu space. 
The effective action S'[p, p^^. A, A^^] is obtained by in- 
tegrating out the fermions. To quadratic order in the 
bosonic fields and their gradient {dr or V), it is suffi- 
cient to retain the first and second order cumulants of S' 



4 



with respect to the mean-field action: 



+ 



dr / d-^r 



2 



5|A 

+ zA?S)|A 



„HS , ■ HS 



8to 2 / 



(16) 



Denoting by and Aq the two eigenvalues of the fiuc- 
tuation matrix appearing in II17|I . we have 



0ci + g 



2£V 



/3q 



(20) 



to order 0{g'^). For g ^ (at fixed a), the mode corre- 
sponding to the eigenvalue Aq is frozen, which leads to 



where the averages (• • • )c are calculated with respect to 
the mean-field action 5mf- Calculating the first and sec- 
ond order cumulants and integrating out the Hubbard- 
Stratonovich fields p^^ and A^^ (Appendix EJ) we ob- 
tain 



dr J d-'rprl^-i 

/3q 



8m 



Mb 
2 



y(dp_„<5|A|_,)( -r-" 



(17) 



where 8pq and (5|A|q are the Fourier transforms of (5pr 
Pv ~ Po and (5|Ar| = |Ar| — Aq, and 



/3q 



ff^7q 



ingo(q)-^ 



.[n-(q)-^n^+(q)]2, 



— n-(q)-ing+(q), 



-n+-(q) - n++(q) + 2n-(q)-ing+(q)'. 



(18) 



We use the notation q — (q, iiOi,) and = w„ where 
is a bosonic Matsubara frequency. The mean-field cor- 
relation function IIqq (q) — {jQ{q)jQ (—<?)) c is calculated 
in Appendix El and n^o"'(q) = ngo'(q,^- = 0). j^iq) is 
the Fourier transformed field of jgr [Eq- CSl]- Eq. ifTzjl 
shows that half the fermion density is the conjugate vari- 
able of the phase 9r of the pairing field. Eqs. II17I18|I 
agree with Eq. (2.3) of Ref. [l3 except for the coefficient 
of 6p-qS\A\q which is found to have opposite sigui^ 

We now discuss the strong-coupling limit (not consid- 
ered in Ref. To leading order in poa"^ and |q|a, we 
have (Appendix |Bt^ 



1 



/3q 
7q 



Airpoa^ 
1 



1/2 



1 + ^TrpoO 



l + 47rpoa' + T|q|'a' 
2pQnia'^ \ 4 

m / 3 ^ 7 , ,o . 



(19) 



6pr 



1 



g \pQm^a 



1/2 



(21) 



Density (Spr) and modulus ((5|Ar|) fiuctuations do not 
fiuctuate independently in the low-energy limit but are 
tied by the relation (|2iJ. From H17I19I21II . we deduce 
that the dynamics of the Fermi superfiuid is determined 
by the effective action 



S[p, 6] = / dr d^r 



T^a 2 



PA - 

{V5p,f 
32pom 



8m 



Pb 
2 



(22) 



Introducing the bosonic field 



(23) 



we recover the standard action of a Bose superfiuid, 



dr / d'^r 



i^ri 9t - PB - 1 Ipr 



2m f 



2iTaB 

TUB 



(24) 



where = 2m and as = 2a are the mass and the 
scattering length of the bosons. The result as = 2a 
corresponds to the Born approximation for the boson- 
boson scattering, while the exact result is as = 0.6a,— 
Eqs. ^ and (I22II are equivalent in the hydrodynamic 
regime where (Vpr)^/po — (Vpr)^/pr'^ 

Thus, we have shown how, by introducing the physical 
fields pr and Ar from the outset and expanding about 
the mean-field state in the strong-coupHng limit, one ob- 
tains the standard action of a Bose superfiuid. Our ap- 
proach should be contrasted with a number of previous 
works^ii^ii^ii^iiSiiii^ where only the pairing Hubbard- 
Stratonovich field A^^ is considered and the expansion 
is carried out about the non-interacting state, which gives 



the action but for the field ^Pr/2e' 
V' field defined in 



instead of the 
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III. LATTICE MODEL 

In this section, we consider the attractive Hubbard 
model on a bipartite lattice, with Hamiltonian 



H = (4cr'+h.c. 



^clcr-U22n-r]n-ri. (25) 



The operator cj^ (cru) creates (annihilates) a fermion 
with spin a at the lattice site r, Cr = (cr|,Crx)"^, and 
rira = cl^Cra- (f; f') denotes nearest-neighbor sites. The 
chemical potential /i fixes the average density po (i.e. the 
average number of fermions per site) and — C/ (C/ > 0) is 
the on-site attractive interaction. 

We are interested in the strong-coupling limit U ^ 
t where fermions form tightly bound composite bosons 
which behave as local pairs. The latter Bose condense at 
low temperature giving rise to superfluidity. In order to 
derive the low-energy effective action, we could follow the 
procedure used in Sec. ^ Here, we shall use a different 
method, based on the mapping of the attractive Hubbard 
model in the strong-coupling limit onto the Heisenberg 
model in a uniform magnetic fieldfSfi, Thus this approach 
is based on a t/U expansion about the t = limit rather 
than on an expansion about the mean-fleld statei^ 

Under the canonical particle-hole transformation^^ 



(-ir 



(26) 



the Hamiltonian becomes (omitting a constant term) 



(r,r') 



clcr' + h.C.) — y^^ t 



+ hod' ] c, 



(27) 



and corresponds now to the repulsive Hubbard model in 
a magnetic fleld hp = h^z along the z axis. 



U 



(28) 



coupled to the fermion spins. The chemical potential 
U/2 in lf27jl . together with particle-hole symmetry, im- 
plies that the system is half-fllled. The density and pair- 
ing operators transform into the three components of the 
spin density operator: 



(29) 



The equation fixing /i, (cjcr) = po, becomes an equation 
fixing the magnetic field: (cJcr^Cr) = po — 1- 

In the strong-coupling limit U ^ t, the Hamiltonian 
l|27ll simplifies into^ft 



H = J E Sr • Sr' - 2ho • ^Sr, 
(r,r') r 



(30) 



where J = At^ /U and Sr is a spin-^ operator. Using 
spin-i coherent states |J7r) (J^r = l)^the action of the 
Heisenberg model l(3(Hl can be written as 



dr 



[(J7r|S^r)-ho-f2r] + J 



{r,r'> 



(31) 

where \ flr) = dr\^r)- 

The effective action S[p, A] of the superfiuid system 
is obtained by rewriting the action ll3T|l in terms of the 
density and pairing fields of the attractive model. In the 
strong-coupling limit, Eqs. (I29II (written now for fields 
rather than operators) become^S, 

Pr = f^r + l, 
A, = — 



2 ' 



(32) 



where — ± zfi^. The condition = 1 implies 
that Pr and Ar do not fluctuate independently but are 
tied by the relation 



|ArH^[Pr(2-Pr)]^/'. 



(33) 



In the low-density limit (pr <C 1), where the Pauli 
principle (which prevents two composite bosons to oc- 
cupy the same site) should not matter, we expect to re- 
cover the stand ard a ction of a Bose superfluid. In that 
limit, |Ar| ~ -y/ Pr/2; the pair density |Arp equals half 
the fermion density pr, and Ar = |Ar|e'^"' coincides with 
the bosonic field ipr — \J pr/2e*^''. To order 0{pf), we 
deduce from H31l.'^3ll 



S[P,0] = 



J 



E 



PrPr 



)] • (34) 



(r,r') 

-(PrPr')'^'(2-Pr)' 



The term {i/2)pr9r comes from the Berry phase term 
(Jlrl^^r) of the action S[fl\ [Eq. Il31|l ] with a proper gauge 
choice. If we further assume that pr and 9r are slowly 
varying in space, we obtain 



S[p,9] = / dr\j2 



i ■ /, Jz\ Jz 2 

-Prfc'r - I fto + — IPr + — Pr 



E(^■■Pr')'/'cOs(0r-er')|• 
(r,r'> ^ 



(35) 
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or, equivalently, 



-tB X! (CV'r' +C.C.) L 
{r,r') ^ 



(36) 



where is — J/2, Ub — 2Jz, = 2/io + and z is 

the number of nearest-neighbor sites. We therefore ob- 
tain the action of the Bose-Hubbard model with on-site 
repulsive interaction Ub and nearest-neighbor hopping 
amplitude ts- In the continuum limit and for a cubic 
lattice, the latter gives a boson mass tub — 1/J as ob- 
tained in Ref. I20- 



IV. CONCLUDING REMARKS 

In this paper, we have shown that a Fermi superfluid in 
the strong-coupHng limit, where superfluidity originates 
from BEC of composit e bos ons, can be described by the 
complex field ipr = \J Pr/2e*^'', where pr is the fermion 
density and 6'r the phase of the pairing field Ar. Such 



is description is made possible by the fact that density 
(pr) and amplitude (|Ar|) fluctuations are not indepen- 
dent in the strong-coupling Hmit. The effective action 
S\p, &\ is derived by introducing the physical fields pr 
and Ar from the outset by means of Lagrange multi- 
plier fields and A^^. The latter play the role of 
the Hubbard-Stratonovich fields usually introduced via 
a Hubbard-Stratonovich transformation of the fermion- 
fermion interaction. 

For continuum models, the effective action is derived 
from an expansion about the mean-field state. It corre- 
sponds to the usual action of a Bose superfiuid of den- 
sity Pr/2 where the bosons have a mass ttib = 2m and 
interact via a contact potential with amplitude qb = 
AiraB/mB, as = 2a. 

For lattice (Hubbard) models, the effective action is ob- 
tained from an expansion about the t = Hmit, using the 
mapping of the attractive Hubbard model in the strong- 
coupling limit onto the Heisenberg model in a uniform 
magnetic field. The effective model is a Bose-Hubbard 
model with an on-site repulsion Ub = 2Jz (with z the 
number of nearest-neighbor sites) and a nearest-neighbor 
intersite hopping amplitude Ib = J/2, where J = At^ /U . 



J 



APPENDIX A: LOW-ENERGY EFFECTIVE ACTION S[p, A] 

In this appendix, we derive the effective action ifTTIl starting from Eq. ifTfiji . The first and second order cumulants 
are given by 



-P^^ n-(,)pf - 2p« ng+(,)Af - 2pH n-(g)AHf - a« n++(<z)A 



-Af X-o-('?)A«r - 2AHSn+-(,)AHf 



(Al) 



The second order cumulant is written in Fourier space. The mean-field correlation function Hqq (q) — (jp (9)jo (~9))c 
is defined in Sec. HlBl To obtain iPTjl . we have used the fact that mean-field correlation functions involving the 
current vanish (Appendix E|- In the low-energy limit, we can approximate Hqq (q) by its static limit Hqq (q) = 

n[;o-'(q,^. = o). 

Integrating out the Hubbard-Stratonovich field p^^, we obtain 



dr / d^r 



-g6\A.f-lSpl 



E{|n„V(q)(A«^Af -f c.c.) -f Af *n„V(q)A;; 



8m 2 ) 

HS 



i{A 



A^^)S\A, 



+ 2noo(q) '[Sp-,Sp,-2iTloUl)Sp-giAf +A_g ;--oovm; V--,-, 



,HS*^ 

-q I 



n;^o+(q)^(A"^A 



HSaHS I aHS*aHS 

-9 9 



2iAf n 



■(A2) 



where 5pr = Pr — Po and (5|Ar| = Ar — Aq. Here we have neglected constant terms and use the saddle-point equations 
To obtain the action S[p, A] in terms of the physical fields only, one has then to integrate out the Hubbard- 
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Stratonovich field A^^: 

[ I?[A«S] exp{-5o[A«S] - ^ ^(Af + A^f )[<5|A|_, - ngg(q)-^ng+(q)5p_j} 

= ^^p{-^ E«^f + ^''r)(^"' + Af *))o[5|A|, - n-(q)-in^+(q)5p,][<5|A|_, - n-(q)-in^+(q)5p_,]} 

= exp{i ^[Mn(q) + + 2MM][S\AU ~ n-(q)-^n^+(q)5p,][5|A|_, - n-(q)-^n^+(q)<5p_,]}, (A3) 

where averages (• • • )o are taken with the Gaussian action 




In the following, we denote by and Bq the diagonal and ofl-diagonal components of M ^(q), and Cq = + Bq. 
The efifective action S[p,A] deduced from Eqs. ljA2IA3|l is gi ven by Eg. (( TTjl. 



APPENDIX B: MEAN-FIELD CORRELATION 
FUNCTION 



In this appendix, we calculate the mean-field correla- 
tion function W^^^riq) = {j^iq)j^'{-q))c {v,v' = x,y,z- 
= 0,x,y,z) in the strong-coupling limit poa^ <^ 1 
and for |q|a ^ 1. 



1. General expression 



j^{q) is the Fourier transformed field oi j^^ [Eq. lfT5|l ]: 



fermionic Matsubara frequency. We have 



~^T.^G{k)Gik + q)-Fik)Fik 



^^G(fc)G(-fc-q), 
I3V ^ m2 {"^ ^ 2 



M 2 



x[G{k + q)G{k) + F{k + q)F{k)], (B2) 

where G and F are the mean-field propagators [Eq. Q]. 
The correlation function {j'^{q)jo{—q)) (m 0) vanishes. 

In the following, we consider the static limit IIqq (q) — 
IIqq (q, uJi, = 0) . Performing the sum over Matsubara 
frequency in ljB2ll in the T = limit, we obtain 



ngo^(q) = 



noo+(q) 



noV(q) 



1 



k -E-k + -Ek+q 

noo"(q) 



^ 0131) n++(q) = noT(q) = -I 



CkCi 



k+q 



£^k-Bk+q 

A^Ck 



£^k-Bk-t 



k (-^k + £^k+q)£'ki?k-(-q 



AHS^ 



k 2(£'k + -Ek+q)-E'k-Ek+q ' 
, CkCk+q 



k2(£;k + £;k+q) 



i^k-Bk+q 



(B3) 



where V is the volume of the system and V ^ J2k ^ /k 
for V oo. A; — (k, iw) and X]fc = X^kw where is a 



The correlation function 11°^, (q) vanishes for q = 0. 
Since j^^ multiplies 9^0r in the action S", it is sufficient 
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to consider H^^, (q = 0) to obtain the eflective action 
S[p,A] to order (a^6'r)^ 

We next expand the correlations to order Odqp). 
Writing ^k+q = Ck + ^k,q with Xk,q = k-q/TO+|qp/2m, 
we obtain 



ngo(q) 



n^o^(q) 



noV(q) 



AT 



3a 



2^5 l^k.q 



2El El 2£;7y k,q 



-L 



ATCk 

k 

3El 



3gk 



k,q 



X, 



k,q 



n++(q) = noo-(q) 



AT 



i-|^^k,q 



1 

2jk 

+ 



4£;^ 
1 

2E^ 
1 

8^ 



k v2 ^ 
^k,q^ 


) 




/a 








sen 


8i?r 





3e 



iEl 



^k,q 



2. Strong-coupling limit poa^ 1 



(B4) 



Other useful relations are obtained by differentiating 
Eqs. ljB5|l with respect to /i. Note that A is sent to in- 
finity whenever the integral over k converges. In the 
strong-coupling limit, the small parameter expansion is 
~ Poo . Approximate expressions of the 
mean-field correlation functions can be obtained by ex- 
panding Eqs. lO and ljB4p in power of AHS^ and usmg 
Eqs. ljB5|l (as well as those obtained from ljB5ll by dif- 
ferentiating with respect to fi). A straightforward (but 
somewhat lengthly) calculation then gives Eq. ^ and 



ng^(q) 



ngo^(q) 



noV(q) 



2 / 1 A 3 |q|^a^ 

poma 1 — 47rpoa 



2 \ 1/2 

Pqiti a ^ 
An 



, 7 3 |q|'a2 

1 Trpna - — 

4 ^' 12 



poma^ . 3 5 I ,2 2 
^(1 -47rpoa - Yg|q| a 



nJo"(q) = - 



1 pQiria 
9 ~ 



ma q ■ 

(1 - inpoa') - —P- 



(B6) 



In the strong-coupling limit, the chemical potential /i 
is negative. We then have 



1 toA TO^/2|/i|l/2 

k^ " ^ \/27r ^ 
dk 



€ 2x/2|m|i/2- 



(B5) 
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